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The masses and dispersions of light hadrons are calculated in lattice QCD 
using an 0(a 2 ) tadpole- improved gluon action and an 0(a 2 ) tadpole-improved 
next-nearest-neighbor fermion action originally proposed by Hamber and Wu. 
Two lattices of constant volume with lattice spacings of approximately 0.40 
Qh, fm and 0.24 fm are considered. The results reveal some scaling violations at 

(— 1 \ the coarser lattice spacing on the order of 5%. At the finer lattice spacing, 

the calculated mass ratios reproduce state-of-the-art results using unimproved 
actions. Good dispersion and rotational invariance up to momenta of p a ~ 1 
are also found. The relative merit of alternative choices for improvement 
operators is assessed through close comparisons with other plaquette-based 
tadpole-improved actions. 
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I. INTRODUCTION 



Lattice discretization of the continuum QCD action introduces errors at finite lattice 
spacing a. The standard Wilson gauge action has 0(a 2 ) errors and the standard Wilson 
fermion action has 0(a) errors. Simulations using these actions have shown that lattice 
spacings of 0.1 fm or less and lattice volumes of 24 4 or larger are needed in order to hold 
systematic errors to the 10% level. Such simulations are major undertakings and require 
enormous computing power to extract even the most basic of hadronic observables, the 
hadron masses (see [|J for an example and for a recent review). 

During the past few years, considerable efforts have been devoted to improving lattice 
actions. The idea is to reduce or remove the discretization errors from the actions so that 
they have better continuum- like behavior. At the same time, errors due to the lattice 
regularization are accounted for through the renormalization of the coefficients multiplying 
the improvement operators. The hope is the use of improved actions will allow one to 
simulate efficiently and accurately on coarse lattices, such that computer resources may be 
redirected to the simulation of QCD rather than quenched QCD. Moreover, one may turn 
the focus of investigation towards quantities of experimental interest. 

The way to develop improved discretizations of continuum actions is not unique. There 
are numerous programs in this endeavor having their foundation in one of two formulations. 
One is based on lattice perturbation theory which is used to derive the renormalized coef- 
ficients through a mean-field reordering of the perturbative expansion, as recently reviewed 
by Lepage 0. The other is based on renormalization group theory to construct fixed-point 
(or perfect) actions, as recently reviewed by Hasenfratz ||. 

In the pure gauge sector, the 0(a 2 ) tadpole-improved action of |5j leads to dramatic 
improvement in the static potential and glueball masses up to lattice spacings of 0.4 
fm. Excellent scaling and topological properties have been demonstrated in fixed-point 
actions 0. In the light quark sector, hadron spectroscopy has been investigated with a 
variety of improved actions including the 0(a)-improved SW action ||, the 0(a 2 )-improved 
D234 action |9j and its variants [p],|10| , and the D%34 action of Hamber and Wu [p4] -|T4]| 



considered here. Fixed-point fermion actions have been studied in free field theory [|T5|J16 
and in simulations 



In the present work, we report further calculations of the light hadron spectrum using 
the 0(a 2 ) tadpole-improved gluon action of || and the D%34 action proposed by Hamber 
and Wu [|11[] many years ago. This action is an 0(a 2 ) next-nearest-neighbor fermion action 
with tadpole- improved estimates of the coupling renormalizations. This action is selected 
primarily due to its simplicity. The cost of simulating it is about a factor of two as compared 
to standard Wilson fermions. Our goal is to study its feasibility as an alternative action to 
SW which has the clover term, or to D234 which has both next-nearest-neighbor couplings 
and the clover term. 



This work extends previous studies of the D%34 action fT3] , |T4|| through finer lattice spac- 
ings, improved statistics, and the simulation of additional observables. In particular, we 
examine dispersion relations and test the rotational symmetry of both the gauge and the 
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fermion actions. Hadron mass ratios are calculated for a wide variety of hadrons including 
hyperons. To explore scaling violations, we consider two coarse lattices of approximately 
fixed physical volume: 6 3 x 12 at a lattice spacing of 0.40 fm and 10 3 x 16 at 0.24 fm. 



II. IMPROVED LATTICE ACTIONS 



A. The Gauge Action 



The improved gauge action employed in this investigation is given by ||: 

Sg = PY, ^ Tr(l - Ufi - ^ £ ifo, Tr(l - U n ) . (1) 
pi 6 2UM ° rt 6 



The second term removes the 0(a 2 ) errors at tree level. Perturbative corrections are esti- 
mated to be of the order of 2-3% ||. Here, U T ^ denotes the rectangular 1x2 plaquettes. uq 
is the tadpole factor that largely corrects for the large quantum renormalization of the links 
U^,(x) = exp(ig J^ +a ^ A(y) ■ dy). In this calculation we use the mean plaquette to estimate 
Wo, 

n \ 1/4 

«o = [~Re Tr(C/ pl )J , (2) 



and will focus our evaluation of lattice action improvement on other plaquette-based im- 
proved actions. Uq is determined self-consistently in the simulation. An alternative choice 
is to use the mean link in Landau gauge u = (^TtU^lg 
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B. The D%34 Fermion Action 



The improvement program of Sheikoleslami and Wohlert (SW) JTJ] provides a systematic 
approach to the improvement of lattice fermion actions. However, the on-shell improvement 
program leaves some freedom in the relative values of the coefficients of the improvement 
operators. In this investigation, we consider a specific case of the general class of D234 
actions |TD[ in which the improvement parameters are tuned to remove the second-order 
chiral-symmetry-breaking Wilson term. This fermion action may be written 

M Dx34 = m g + 7- V + ^]T (-a 2 V„ A„ + b a 3 Aj) (3) 



where 



V M ip(x) = — U^x) i)(x + fx) - Ul(x - ji) i)(x - /i) 



(4) 



3 



and 



a 1 



Ufj,(x) ip(x + fi) + U^x - fi) ip(x — n) — 2ip(x) 



The second-order term of the D234 action 



a 



(5) 



(6) 



breaks chiral symmetry and does not appear in the D%34 action. However, the fourth-order 
term of (|3]) also breaks chiral symmetry and provides for the removal of the fermion doublers. 

Explicit evaluation of (|f) combined with a Wilson-fermion-style field renormalization 
factor of (3 k/2) discloses the following simple fermion action |D]|L2"fl : 

S F = J2^( x )^( x ) 



U (r) - UUx) 

${x) (b - 7a4 ) ^{x + ») + $(x + m) (b + 7 „) ^x) 



u 



u 



Wo 



-H>{x + 2/x -- b - - 7m ^ il>{x) 

4 8 Uq Uq 



(7) 



The tadpole-improvement factors are explicit here. These coefficients remove both 0(a) and 
0(a 2 ) errors at tree level. The bare quark mass is related to k and b through 

2 



3k 



-46. 



Thus the renormalized quark mass is given by 

2 n 



m q ,-, . 

The wave function renormalization factor is also modified: 

/3 K 



cont. 




latt. 



(8) 



(9) 



(10) 



The standard Wilson action can be recovered by truncating the two-link terms of the action, 
redefining k and renormalizing the fermion field operators. 

Free dispersion relations can be obtained by locating the poles in the fermion propagator. 
Fig. [I] shows free dispersion relations for the continuum, Wilson, and D%34 fermions. It is 
clear that D%34 fermions follow the continuum more closely than Wilson fermions. Note 
that there exists an unphysical high energy doubler (or ghost) in the D%34 action. It is very 
similar to that for the D234 action || and is a general feature of fermions with next-nearest- 
neighbor couplings. The doubler can be 'pushed away' from the low momentum region by 
various techniques, such as tuning the value of b, or using an anisotropic lattice |1(J. We 
simulate with 6=1, which gives good dispersion to pa ~ 1 as illustrated in Fig. |T[ 
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0.0 



pa 



FIG. 1. Free dispersion relations for Dx34, Wilson, and continuum fermions. Momentum p 
is along the (1,1,0) direction. Beyond the D%34 branch point, the real part of the two conjugate 
roots is shown. The upper graph is for bare mass ma = 0, the lower ma = 0.3. 
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III. LATTICE SIMULATIONS 



A. Methods and Parameters 



Quenched gauge configurations were generated using the Cabibbo-Marinari pseudo-heat- 
bath method Periodic boundary conditions were used in all directions for the gauge 



field and in spatial directions for the fermion field. Dirichlet boundary conditions were used 
for the fermion field in the time direction. Configurations separated by 300 sweeps were 
selected after 4000 thermalization sweeps from a cold start. 

Wilson loops were computed using a method developed in Ref. ||20|| . Both on-axis and off- 
axis loops were considered. The static potential was determined by fitting exp(— V(R) T) = 
W(R, T). The string tension a was extracted from the ansatz V(R) = Vq + crR — E/R where 
Vq and E are constants. Fig. ^ shows our results. Good rotational invariance of the static 
potential is observed. Table | gives the extracted string tension and the Rq parameter, 
defined by R 2 dV(R)/dR = 1.65 at Rq [F2T]]. Its physical value is R ~ 0.50 fm and it is 
independent of the quenched approximation. Using ^fa = 440 MeV to set the scale, the 
lattice spacings were determined as 0.40(6) fm and 0.23(1) fm, respectively. Using Rq to set 
the scale, they are 0.40(3) fm and 0.22(1) fm. Rq depends less on the functional form of 
V(R) than the string tension. 





TABLE I. String 


tension and Rq extracted from 


the static potential. 




p 


aa 2 


Ro/a R m i n /a 


Rmax 1 Q> 


X*/Ntsf 


6.25 


0.81(23) 


1.26(11) 1.00 


4.24 


1.46 


7.00 


0.254(8) 


2.279(18) 1.73 


6.93 


0.91 



Five quark propagators were computed by the Stabilized Biconjugate Gradient algo- 
rithm []22] for each configuration. The five k values correspond to quark masses of roughly 
210, 180, 150, 120, 90 MeV, respectively, for both lattices. The second value, 180 MeV, 
was taken as the strange quark mass. A point source was used at space-time location 
(x,y,z,t) = (l,l,l,2) on the 6 3 x 12 lattice and (1,1,1,3) on the 10 3 x 16 lattice. The gauge- 
invariant smearing method of [23] was applied at the sink to increase the overlap of the 
interpolating operators with the ground states. Table || shows a summary of the parame- 
ters. 

Statistical errors were estimated in a third-order, single-elimination jackknife, with bias 



TABLE II. Summary of lattice parameters for the spectrum calculations. Nu is the number 
of configurations considered. The second kappa was selected as the one that corresponds to the 
strange quark mass. 



Lattice 


P 


Nu 


u 






K 








6 3 x 12 


6.25 


155 


0.820 


0.175, 


0.177, 


0.180, 


0.183, 


0.187 


0.1967(2) 


10 3 x 16 


7.00 


100 


0.866 


0.170, 


0.172, 


0.174, 


0.176, 


0.178 


0.1823(3) 



6 



5.5 



4.5 - 



3.5 



td 



2.5 



1.5 



0.5 




5 

R/a 



10 



FIG. 2. Static potential from Wilson loops. The empty circles are for (3 = 6.25, the solid circles 
for = 7.0. The lines are best fits. The statistical errors were from 200 configurations in both 
cases. 



corrections p4j. A third-order jackknife provides uncertainty estimates for the correlation 
functions, fits to the correlation functions, and quantities extrapolated to the chiral limit. 



B. Effective Masses 



Masses were extracted from lattice correlation functions in several time slice intervals. 
The effective hadron mass as a function of time is defined by 



M(t + 0.5) = log[G(t)] - log[G(t + 1)] 



(11) 



where G(t) is the standard two-point function in Euclidean space-time. As usual, effective 
mass plots for the pion are very flat. All fits of time slices 4 through 10 on our coarse lattice 
and 6 through 12 on our fine lattice agree within one standard deviation. 

Fig. |] displays effective mass plots for the p meson and nucleon. The p meson displays 
good plateau behavior. We find xV-^df ~ 0.5 for correlation function fits from t = 4 
through 10 on our coarse lattice and xV^df — 0.7 for correlation function fits from t = 6 
through 11 on our fine lattice. Again, all fits agree within one standard deviation. 

On our coarse lattice the nucleon mass displays good plateau behavior. The xV^df 
ranges from 1.0 to 0.1 for the lightest quark mass. All fits in the interval 4 through 10 agree 
within one standard deviation. 
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FIG. 3. Effective mass plots for the p and nucleoli. The first column is for the coarser lattice, 
the second for the finer lattice. The five quark masses decrease by value from top down. 

On our fine lattice, plateau behavior for the nucleon mass is excellent within the regime 
t = 6 through 9. For the heaviest three quark masses, however, the xV^df identifies two 
separate fitting regimes providing acceptable fits, namely 6 through 9 and 9 through 12. The 
latter regime disappears for the lightest two quark masses while the former regime grows, 
providing acceptable fits for t — 6 through 10. Since the effective masses rise from below as 
a function of t for smeared sinks, we report results for the regime t — 6 through 9 where the 
signal is less likely to be dominated by noise. 

We note that this drift in the effective mass plot is only significant as measured by the 
X 2 /Ndf for the ground state nucleon mass. We find that other quantities such as the nucleon 
energy or p-meson mass are sufficiently correlated to provide results insensitive to the fit 
regime. For example, dispersion relations agree within one standard deviation for all fits of 
t — 6 through 12 with no systematic drift in the central values. Moreover the extrapolated 
nucleon to rho mass ratios for the time slice regimes 6 through 9 and 9 through 12 are 
1.44(5) and 1.41(13) respectively. Again, the deviation of central values of the distributions 
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is small relative to the statistical uncertainties. 

Effective mass plots for A are acceptably flat but suffer a loss of signal earlier in Euclidean 
time as expected. On our coarse lattice we find acceptable xV^df — 1.0 — * 0.4 for intervals 
t = 4 or 5 through 10 with smaller values corresponding to lighter masses. On our fine 
lattice we find acceptable xV^df — 1-0 — > 0.6 for intervals t = 6 or 7 through 12, again 
with smaller values corresponding to lighter masses. 

In summary, we report values taken from covariance matrix fits to the time slice interval 
4 through 8 on our coarse lattice and 6 through 9 on our fine lattice. These regimes provide 
the best signal-to-noise and good correlated xV-^df- 



C. Hadron Masses 



Fig. [| shows the extracted hadron masses plotted as a function of the quark mass. For 



future reference, these results are tabulated in Table |TTT| . n cr is determined by linearly 
extrapolating m 2 , as a function of m q to zero. The form M = Cq + c\ m q was used for other 
extrapolations to the chiral limit. Fits with an additional term c^vn^j 2 were also considered 
and similar results were found with slightly larger error bars. The correct ordering of all the 
states is clearly seen at both values of (3. 

Ratios of the chirally extrapolated masses are given in Table along with the ratios as 
observed in nature [^J. At (3 — 7.0 the lattice spacing estimates follow the familiar pattern 
having the value based on the string tension lying between that of the p and nucleon based 
values. This is most likely an artifact of the quenched approximation. However, at f3 — 6.25 
we find significant disagreement among the values and an unusual reordering of values. 

Focusing first on ratios of hadrons having the same angular momentum, we see very 
little change in the values as the lattice spacing is decreased. These ratios are also remark- 
ably similar to those observed in nature, despite the fact that these are quenched QCD 
calculations. In addition, these ratios support our selection for the strange quark mass. 

This close resemblance to nature is not shared by ratios of hadrons with different angular 
momentum. All four classes of ratios significantly disagree with those of nature. Once again 
we see the familiar quenched artifact of the Octet/Vector ratio being too large and the 
Decuplet/Octet ratio being too small. 

The standard failure of the K/p mass ratio in the quenched approximation is also 
seen here. This shortcoming has been widely realized through an examination of the J- 



parameter |56| defined by 



J = jji — 

p dml 



dm n rriK* — m p 

m p /m n =i.8 ^ m 2 — rn? 



™K* Z 2 Z± ■ (12) 

K 



Empirically this ratio is 0.48. However we find 0.40 on our coarse lattice and 0.38 on our 
fine lattice. The physics associated with this discrepancy was first reported in Ref. p? 



There it was pointed out that the self-energy generated by two-pion intermediate states 
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-0.1 0.0 0.1 0.2 0.3 0.4 0.5 

m a 

g 




-0.05 0.00 0.05 0.10 0.15 0.20 0.25 0.30 

m a 
i 



FIG. 4. Hadron masses in lattice units as a function of m q . The upper graph is for the coarser 
lattice, the lower for the finer lattice. The lines are chiral fits, linear in the quark mass. For better 
viewing, the decuplet masses are shifted upward by 0.5 unit. 
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TABLE III. Hadron masses as a function of m q (both in lattice units). The results were 
extracted from time slices 4 to 8 on the coarser lattice, and 6 to 9 on the finer lattice. The last 
column reports results after extrapolations to the chiral limit. 

6 3 x 12 lattice at (3=6.25 



Hadron m q =0A2 m q =0.37 771^=0.31 m g =0.25 m q =0.17 m q =0 



7T 


1.321(6) 


1.248(6) 


1.136(6) 


1.020(7) 


0.855(8) 





K 


1.29(1) 


1.25(1) 


1.19(1) 


1.14(1) 


1.07(1) 


0.91(1) 


P 


1.66(1) 


1.61(1) 


1.53(2) 


1.46(2) 


1.36(3) 


1.15(4) 


K* 


1.64(1) 


1.61(1) 


1.57(2) 


1.54(2) 


1.49(2) 


1.39(2) 


N 


2.67(3) 


2.57(3) 


2.43(3) 


2.29(4) 


2.10(5) 


1.70(6) 


A 


2.64(3) 


2.57(3) 


2.48(3) 


2.38(3) 


2.25(4) 


1.98(5) 


E 


2.64(3) 


2.57(3) 


2.48(3) 


2.39(3) 


2.28(4) 


2.02(5) 




2.61(3) 


2.57(3) 


2.52(3) 


2.48(3) 


2.41(3) 


2.27(4) 


A 


2.86(4) 


2.78(4) 


2.66(4) 


2.54(4) 


2.40(6) 


2.07(7) 


S* 


2.83(4) 


2.78(4) 


2.70(4) 


2.62(4) 


2.52(5) 


2.30(6) 


"* 


2.82(4) 


2.78(4) 


2.74(4) 


2.70(4) 


2.65(4) 


2.54(5) 


10 3 x 16 lattice at 0=7.0 












Hadron 


m q =0.26 


m q =0.2l 


m q =0.17 


m g =0.13 


m q =0.08 


m q =0 


7T 


0.967(7) 


0.871(7) 


0.772(7) 


0.667(7) 


0.548(8) 





K 


0.920(7) 


0.871(7) 


0.823(7) 


0.775(7) 


0.726(7) 


0.637(5) 


P 


1.19(1) 


1.12(1) 


1.05(1) 


0.99(2) 


0.92(2) 


0.80(2) 


K* 


1.15(1) 


1.12(1) 


1.08(1) 


1.05(1) 


1.02(2) 


0.96(2) 


N 


1.86(2) 


1.74(2) 


1.62(2) 


1.50(2) 


1.37(3) 


1.15(3) 


A 


1.82(2) 


1.74(2) 


1.65(2) 


1.57(2) 


1.48(2) 


1.33(2) 


£ 


1.81(1) 


1.74(2) 


1.66(2) 


1.59(2) 


1.52(2) 


1.38(3) 




1.78(2) 


1.74(2) 


1.69(2) 


1.65(2) 


1.61(2) 


1.53(2) 


A 


2.00(2) 


1.90(2) 


1.80(2) 


1.71(3) 


1.62(4) 


1.44(4) 


S* 


1.97(2) 


1.90(2) 


1.84(2) 


1.78(3) 


1.72(3) 


1.60(3) 




1.94(2) 


1.90(2) 


1.87(2) 


1.84(2) 


1.81(3) 


1.75(3) 
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TABLE IV. Mass ratios after extrapolation to the chiral limit. a p and a at are lattice spacings 
set by the rho mass (770 MeV) and the nucleon mass (938 MeV). 





jj — O.ZO 


p — / .u 


rL/XpT;. 


a s t 


U.4U(oJ tm 


U.2zU(2j tm 




a p 


0.30(1) tm 


0.205(6) tm 




a N 


0.36(1) fm 


0.242(6) fm 




Vector /Vector 








K*/p 


1.20(2) 


1.20(2) 


1.16 


i 1 

91 P 


1.40(3) 


1.39(3) 


i on 

1.32 


Octet/Octet 








A / AT 
1\/ iv 


l.lO^z J 


1.16(1] 


i.iy 


E/N 


1.19(2) 


1.20(1) 


1.27 


E/N 


1.34(3) 


1.33(2) 


1.40 


Decuplet /Decuplet 








Zj /A 


i io/i\ 

1.12(1) 


1.11(1) 


1.12 


a /A 


1.24(3) 


1.22(2) 


1.24 




1.34(3) 


1.33(2) 


1.36 


Pseudoscalar / Vector 








K/p 


0.79(2j 


0.79(2j 


0.64 


W^tCl / V CtlUI 








N/p 


1.48(6) 


1.44(5) 


1.22 


K/p 


1.72(6) 


1.66(5) 


1.45 


X,/p 


1. (0(0 J 


1. I Z[0 ) 


1.00 


a/p 


1.9/(oJ 


1.91(0 J 


1 71 

1.71 


JJctlipicl/ VcCLUI 








A/p 


1.79(7) 


1.77(5) 


1.60 


S*/p 


2.00(6) 


1.97(6) 


1.80 


s*/p 


2.20(7) 


2.16(6) 


1.99 


n/p 


2.40(7) 


2.35(6) 


2.17 


Decup let/Octet 








A/N 


1.21(5) 


1.23(2) 


1.31 


T,*/N 


1.35(5) 


1.37(2) 


1.47 


E*/N 


1.49(4) 


1.50(3) 


1.63 


n/N 


1.63(4) 


1.64(3) 


1.78 
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1.10 




0.70 1 1 1 1 1 1 1 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

m 2 (GeV 2 ) 



FIG. 5. Plot of the p-meson mass as a function of the squared pion mass obtained from our 
finer lattice. a p has been used to set the scale. The dashed line illustrates the standard linear 
extrapolations of m^. and m p . The solid and dot-dash curves include the two-pion self-energy of 



the p meson [27] for dipole dispersion cut-off values of 1 and 2 GeV respectively. The increase in 
the slope at mp/m^ = 1.8 (m^ ~ 0.21) provided by the two-pion self energy is the right order of 
magnitude to restore agreement with the empirical value. 

of the p-meson, which is excluded in the quenched approximation, acts to increase the J 
parameter. Fig. || provides a sketch of how including the two-pion self-energy of the p can 
increase the value of J from 0.38 to 0.46. 



Perhaps the most important information displayed in Table [TV] is that the Octet/Vector 
mass ratios display less than satisfactory scaling for the larger lattice spacing. To further 
examine scaling and make contact with other studies, we focus on the the N/p mass ratio 
which is among the the most revealing of ratios. 

Fig. ^ shows a comparison of the N/p mass ratio versus M p a at the chiral limit. Fig. [7] 
shows the N/p mass ratio as a function of M p a at a fixed n/p mass ratio of 0.7 f29 |. This 



method is free of complications from chiral extrapolations. Both cases clearly show the 
improvement provided by the D^34 action. Indeed, the D%34 action has reproduced the 
state-of-the-art quenched QCD ratios using unimproved actions at coarse lattice spacings of 



0.24 fm. Our results at the chiral limit are compatible with those in | 13[j , but lie slightly 
above the SW and D234 results. The results from various improved actions at fixed ir/p ratio 
show a certain degree of universality over a wide range of lattice spacings. The performance 
of fixed-point actions at very coarse lattice spacings is remarkable. 
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FIG. 6. The N/p mass ratio versus M p a at the chiral limit. Solid symbols denote the standard 
Wilson action. Open symbols denote improved actions including SW ||] (o), D234 || (A), D^34 
H (□), and D X 34 (o) (this work). 
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FIG. 7. The N/p mass ratio versus M p a at a fixed n/p mass ratio of 0.7 for various actions. 
The solid symbols denote the standard actions: Wilson (square and diamond), staggered (circle). 
The open symbols denote improved actions: nonperturbatively-improved SW J28| (x), fixed-point 



action |7| (□), SW (o), D234 (A), and D X 34 (o) (this work). 
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D. Dispersion and Rotational Symmetry 



In addition to mass ratios, hadron states at finite momentum projections pa = n{2ir / L) 
were also calculated. Dispersion was examined by calculating the effective speed of light, 
c eff , defined by c 2 ff = (E 2 (p) — E 2 (0))/p 2 , which is to be compared with 1. 

A comparison with SW and D234 lattice actions is made in Table |V[ The lattice 
spacings are based on charmonium for SW and D234 actions whereas the static quark 
potential is used for our results. As such, the lattice spacings in Table |V| are approximate. 
Lattice volumes range from 2.0 to 2.4 fm. D%34 results are based on simulations at k = 0.183 
for our coarse lattice and k = 0.176 on our fine lattice. The dispersion for the 0(a)-improved 
SW action is very poor relative to the excellent dispersions of the next-nearest-neighbor 
improved Dx34 actions. The D%34 dispersion is excellent even at our coarse lattice spacing. 



TABLE V. Comparison of SW, D234 and D%34 actions for the speed of light squared obtained 
from the dispersion of tt and p mesons at m^/nip ~ 0.7 for pa = (2n/L). 



Hadron 


a (fm) 


SW 


D234 


D X 34 


TT 


0.40 


0.63(2) 


0.95(2) 


0.99(3) 


TT 


0.24 




0.99(4) 


1.04(4) 


P 


0.40 


0.48(4) 


0.93(3) 


0.93(6) 


P 


0.24 




1.00(6) 


0.99(6) 



Rotational symmetry is explored in Table [Vl, which reports results for m q ~ 180 MeV. At 



the coarser lattice spacing, some drift in the central values is seen for the pion and nucleon. 
The drift in the pion is similar to that seen for the D234c action reported in Ref . [Q] . However, 
the drift in dispersion for the (f> meson reported in Ref. is not apparent in our results for 
the D^34 action. However, a similar drift may be hidden in the uncertainties. The D%34 
action has much better rotational symmetry than the SW action ||. Moreover, the D%34 
action provides satisfactory dispersion at our finer lattice spacing and is competitive with 
the D234 action @. 

TABLE VI. Evaluation of dispersion and rotational invariance via the effective speed of light. 
The results are for m q ~ 180 MeV. 



a (fm) 


n 


TT 





N 


n 


0.40 


(1,0,0) 


0.98(2) 


0.91(4) 


1.00(7) 


0.99(12) 




(LLO) 


0.91(4) 


0.91(6) 


0.94(8) 


0.91(7) 




(1,1,1) 


0.86(9) 


0.92(10) 


0.92(7) 


0.90(10) 


0.24 


(1,0,0) 


1.04(3) 


1.02(4) 


1.10(6) 


1.06(7) 




(1,1,0) 


1.05(4) 


1.02(4) 


1.06(5) 


1.11(5) 




(1,1,1) 


0.98(6) 


0.98(6) 


1.06(6) 


1.06(5) 
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IV. CONCLUSION 



We have computed masses and dispersion relations of light hadrons in lattice QCD using 
tree-level 0(a 2 ) tadpole-improved gauge and fermion actions. Compared to standard Wilson 
actions, the gluon action has an additional six-link (rectangle) term while the D%34 fermion 
action has an additional two-link (next-nearest-neighbor) term. These actions have the 
appeal of being simple to implement and inexpensive to simulate. 

A great deal of effort is being directed toward finding the ultimate improved action that 
will facilitate simulations on the coarsest of lattices. We note however, that many quantities 
of phenomenological interest such as hadron form factors involve momenta on the order of 
a GeV. As such, a highly improved action which is costly to simulate may not be the ideal 
action for hadron phenomenology, especially for exploratory purposes. 

The mass ratios obtained from the D%34 action at 0.24 fm on a modest 10 3 x 16 lattice 
reproduce the state-of-the-art results using conventional unimproved actions. Excellent dis- 
persion and rotational invariance up to pa ~ 1 are also found. These results demonstrate 
that the D%34 action can serve as a viable candidate for the study of hadron phenomenology, 
and in our view is preferable to the highly-improved but more costly D234 action. We plan 
to use the D%34 action to study hadron properties beyond the spectrum, such as multipole 
form factors of hadrons in general. These results also bode well for future explorations 
beyond the quenched approximation. 
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